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Abstract 

We report the discovery of a significant excess of candidate Ha emitters (HAEs) in the field of the 
radio galaxy 4C 23.56 at z = 2.483. Using the MOIRCS near- infrared imager on the Subaru Telescope we 
found 11 candidate emission-line galaxies to a flux limit of ^ 7.5 x 10~^^erg s~^ cm~^, which is about 5 
times excess from the expected field counts with ^ 3-cr significance. Three of these are spectroscopically 
confirmed as redshifted Ha at z = 2.49. The distribution of candidate emitters on the sky is tightly confined 
to a 1.2-Mpc-radius area at z = 2.49, locating 4C 23.56 at the western edge of the distribution. Analysis 
of the deep Spitzer MIPS 24 /im imaging shows that there is also an excess of faint MIPS sources. All but 
two of the 11 HAEs are also found in the MIPS data. The inferred star-formation rate (SFR) of the HAEs 
based on the extinction-corrected Ha luminosity (median SFR ^ 100 M© yr~^) is similar to those of HAEs 
in random fields at z ^ 2. On the other hand, the MlPS-based SFR for the HAEs is on average 3.6 times 
larger, suggesting the existence of the star-formation significanly obscured by dust. The comparison of 
the Ha-based star-formation activities of the HAEs in the 4C 23.56 field to those in another proto-cluster 
around PKS 1138-262 at z = 2.16 reveals that the latter tend to have fainter Ha emission despite similar K- 
band magnitudes. This suggests that star- formation may be suppressed in the PKS 1138-262 protocluster 
relative to the 4C 23.56 protocluster. This difference among the HAEs in the two proto-clustcrs at z > 2 
may imply that some massive cluster galaxies are just forming at these epochs with some variation among 
clusters. 

Key words: Galaxies: protocluster — Galaxies: individual (4C 23.56) — Galaxies: high redshifts — 
Galaxies: Ha emitters 



1. Introduction 

Clusters of galaxies are the most massive bound sys- 
tems in the universe. The present-day cluster popula- 
tion is dominated by early-type galaxies which make up a 
tight "red sequence" in color-magnitude diagram (Dressier 
1980). The analyses of the zero-point, slope, and scatter of 
the color-magnitude relation have been playing key roles 
in understanding when and how galaxies in clusters are 
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formed and evolved (e.g.. Bower Lucey, & Ellis 1992). 

Studies of low to intermediate redshift clusters indicate 
that the majority of stars in cluster galaxies are formed 
at z > 2 (e.g., Elhs et al. 1997; Bower et al. 1998; Kodama 
et al. 1998; Stanford et al. 1998; Holden et al. 2004; Tran 
et al. 2007). On the other hand, detailed morphological 
analyses of the cluster red galaxies suggest that a signifi- 
cant fraction of the present-day cluster early-type galaxies 
must have been fully assembled after z 1 (e.g.. Dressier 
et al. 1997; van Dokkum et al. 2000; van Dokkum & Franx 
2001; Ellis et al. 2006; Mei et al. 2009). A significant rise in 
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star-formation activity in the cluster cores is also seen as 
one goes to higher (z > 0.4) redshift, the "Butcher-Oemler 
effect" (Butcher & Ocmler 1978; Butcher & Oemler 1984). 

Studies of the Butcher-Oemler effect, as well as the 
morphological transformation in cluster cores, have been 
among the major topics in extragalactic astronomy for 
the past two decades (e.g.. Dressier et al. 1994; Couch et 
al. 1994; Dressier et al. 1997; Balogh et al. 1997; Oemler 
et al. 1997; Poggianti et al. 1999; van Dokkum et al. 
2000; Ellingson et al. 2001; Kodama & Bower 2001; Mei 
et al. 2009; Pannella et al. 2009). As the observational 
data have been accumulated for ever more distant galaxy 
clusters, the "break" of the red sequence at fainter magni- 
tudes and its evolution has attracted great attention (e.g., 
Kajisawa et al. 2000; Kodama et al. 2004; De Lucia et al. 
2004; Tanaka et al. 2005; De Lucia et al. 2007; Gilbank 
et al. 2008). These phenomena have been discussed in 
the framework of the "down-sizing" of the star-formation 
activity, whereby the stellar mass of the galaxies which 
host star formation at a fixed rate is, on average, decreas- 
ing with cosmic time. The evolution of the faint end of 
the red sequence is also connected to the continuous infall 
of star-forming galaxies into the cluster cores, their sub- 
sequent cessation of star-formation, and change in their 
morphologies (Poggianti et al. 2004; Bundy et al. 2006). 

From the local universe to intermediate redshifts (z < 
1), a number of studies suggest that star formation is 
strongly suppressed in the cluster cores compared with 
the low-density field environment (e.g., Hashimoto et al. 
1998; Gomez et al. 2003; Nakata et al. 2005; Poggianti 
et al. 2006; Patel et al. 2009). For passively-evolving 
early-type galaxies, there is also weak, but growing, evi- 
dence that the luminosity-weighted ages of galaxies tend 
to be older in higher density environments. This suggests 
that perhaps the time scale of star formation and/or sec- 
ondary episodes of the star formation vary with environ- 
ment (Thomas et al. 2005; Clemens et al. 2006; Rakos 
et al. 2008; Clemens et al. 2009; Blanton & Moustakas 
2009; Pannella et al. 2009; Thomas et al. 2010). Recently, 
however, this trend is observed to flatten or even reverse, 
with the density dependence of star-formation disappear- 
ing or with enhanced star-formation in the cluster core at 
z> 1 (Elbaz et al. 2007; Ideue et al. 2009; Hayashi et al. 
2010; Tran et al. 2010). 

If the majority of the massive cluster galaxies {M^, > 
lO^*' Mq) begin forming their stars at z > 2, they will 
have high star formation rates at higher redshifts. Hayashi 
et al. (2010) have surveyed star-forming galaxies in the 
X-ray cluster XMMXCS J2215.9-1738 at z = 1.46 us- 
ing a narrowband filter that samples the redshifted [O II] 
doublet located at the cluster redshift. They showed 
that the [O II] emitting star-forming galaxies are preva- 
lent everywhere in the cluster including its core region 
(see also, Hilton et al. 2010). Although star formation 
is still confined to the cluster outskirts in other z ~ 1.4 
cluster (Lidman et al. 2008; Bauer et al. 2010; StrazzuUo 
et al. 2010), the widespread star formation in XMMXCS 
J2215. 9-1738 is in strong contrast to similar surveys of 
lower redshift clusters at z = 0.81 and z = 0.41 (Kodama 



et al. 2004; Koyama et al. 2010). The color- magnitude 
relation for XMMXCS J2215.9-1738 also shows a larger 
scatter than those of other intermediate redshift clusters 
(Hilton et al. 2009; Papovich et al. 2010). Although we 
need a larger sample of clusters at intermediate and high 
redshifts to discern a general trend in the star formation 
— galaxy density relation (Elbaz et al. 2007; Poggianti et 
al. 2008; Sobral et al. 2010), the existence of widespread 
star formation activity in this z = 1.5 cluster strongly sug- 
gests that we are gradually approaching the major epoch 
of galaxy formation in cluster environments. 

At higher redshifts (z > 2), the number of the known 
(proto-) clusters is still quite small, and the majority 
of the known proto-clusters were discovered in narrow- 
band Lya-emitter surveys. With some exceptions (e.g., 
Steidcl et al. 1998; Steidel et al. 2000; Shimasaku et 
al. 2003; Steidel et al. 2005; Gobat et al. 2010), much 
of the search for proto-clusters has been done around 
powerful radio galaxies and QSOs with known redshifts 
(e.g.. Keel et al. 1999; Pentericci et al. 2000; Kurk et al. 
2000; Yamada et al. 2001; Venemans et al. 2002; Kurk 
et al. 2004a; Venemans et al. 2005; Venemans et al. 
2007; Matsuda et al. 2009). The success of these searches 
around radio galaxies is related to the fact that the host 
galaxies of the powerful radio sources tend to be the most 
massive galaxies at any epoch of the universe (De Brcuck 
et al. 2002; Seymour et al. 2007). Powerful radio sources 
are powered by central super-massive black holes, and the 
large inferred masses of the black holes in turn imply mas- 
sive host stellar systems (Magorrian et al. 1998). The 
largest galaxies in the universe tend to lie in the centers 
of largest-scale haloes which must have originated from 
the primordial density peaks where galaxy formation pro- 
cesses are accelerated in a biased way (Rccs 1985). 

Recently Tanaka et al. (2010) analyzed the broad- 
band SED of galaxies in the PKS 1138 proto-cluster and 
found that proto-cluster galaxies show suppressed star- 
formation on average compared to field galaxies (see also, 
Zirm et al. 2008; Doherty et al. 2010). This implies that 
even at z 2 the activity in the cluster core region seems 
to lag behind the star-formation activity we see in galax- 
ies in field surveys. Steidcl ct al. (2000) also reported 
older stellar ages for their z = 2.3 proto-cluster than for 
field galaxies in the same survey, implying a similar trend. 
The detection of a possible (proto-) red-sequence in four 
known proto-cluster fields at z ^ 2 via JHK photome- 
try (Kodama et al. 2007) also indicates that, though the 
star-formation is also likely to be prevalent there, some 
X-bright galaxies may already be converging to the red- 
sequence. 

At higher redshifts, there are some hints that more in- 
tense star-formation is occurring even in the proto-cluster 
environment. In the SSA22 field, which is known to con- 
tain a large overdensity of galaxies at z = 3.1 (Steidel et 
al. 1998; Steidel et al. 2000), Matsuda ct al. (2004) de- 
tected a significant number of the extended Lya blobs 
which are thought to be progenitors of massive galaxies 
(e.g., Dijkstra & Loeb 2009). Two of the largest Lya blobs 
are identified with red, massive (M^, ~ 10^^ Mq) galax- 
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ies, comparable to dominant cluster galaxies at z < 1.5 
(Uchimoto et al. 2008; Stott et al. 2010). Tamura et al. 
(2009) have detected an excess of AzTEC 1.1 mm-sources 
in the same SSA22 field. The positional correlation of 
the AzTEC sources with the Lya emitter distribution 
at z = 3.1 suggests that dust-enshrouded vigorous star- 
formation might be going on in the field. There is also a 
report of a significant excess of submm-bright objects in 
the TNJ 1338-1942 field at z = 4.1, which is known to have 
a proto-cluster identified by Lya emitters (De Breuck et 
al. 2004; Venemans et al. 2007). 

As was demonstrated by Kodama & Bower (2001), star 
formation activity in cluster cores should have been as 
vigorous as in the field environment at some point in the 
past (and may exceed it at even higher redshifts). To ob- 
serve this crossover point we need larger samples of proto- 
clusters at z =2-4. 

The field around 4C 23.56 is one of the most promis- 
ing proto-cluster candidates at z > 2 (Knopp & Chambers 
1997, hereafter KC97). The authors detected an excess of 
candidate Lya emitters as well as an excess of very red 
objects in I — K' color in their survey field (1.25 arcmin^). 
More recently, Kajisawa et al. (2006a) took deep J and 
Kg images with CISCO on Subaru and found a signifi- 
cant excess of Distant Red Galaxies (DRGs, Franx et al. 
2003) in a 1.6 x 1.6 arcmin^ field. This further supports 
the idea that there are a number of evolved galaxies in 
this proto-cluster at z '--^ 2.49. However, both of these sur- 
veys have very limited field coverage and they may not be 
probing the whole picture of the proto-cluster. So far no 
other follow-up observations, including any redshift con- 
firmation, have been reported in the literature. 

Here we present the result of our new Ha emitter 
(HAE) search around 4C 23.56 with the wide-field near- 
infrared imager, MOIRCS, mounted on the 8.2m Subaru 
Telescope. In this study, we show that there is a clear 
and strong excess of the HAEs near the radio galaxy. We 
present their star formation properties and discuss the 
characteristics of this proto-cluster by comparing with 
other fields at similar redshifts. Archival Spitzer MIPS 
data are also used to probe dusty star forming activity in 
the proto-cluster. This is the first narrow-band Ha emit- 
ter survey of a proto-cluster field at z > 2 after Kurk et al. 
(2004a). The field size we cover is nearly 19 times larger 
than that of KC97, and is more than twice as large as that 
of Kurk et al. (2004a). We note that the survey redshift 
of z ~ 2.5 is close to the red end of the Ha-observable win- 
dow (z < 2.7) from the ground, due to the significant rise 
of the thermal background beyond the K-hand window. 

Throughout the paper, the standard ACDM cosmology 
with fi™ = 0.3, f^A = 0.7, h = 0.7 is assumed (1' = 1.18 
comoving Mpc at z = 2.48). All the magnitudes presented 
in this paper are in the AB system unless otherwise noted. 

2. Observations and Data Reduction 

2.1. Observations 

The primary data for the current work were col- 
lected using the Multi-Object Infrared Camera and 



Spectrograph (MOIRCS, Ichikawa et al. 2006; Suzuki et 
al. 2008) mounted on the Cassegrain focus of the Subaru 
8.2-m Telescope at Mauna Kea. The observations were 
executed during the second half-night of UT 1st June, 
2007, during the author's (IT) Subaru Staff Time. The 
weather was photometric but with variable seeing in the 
range 0."47-0"7. 

MOIRCS has two independent sets of optics (including 
both filters and detectors), which divide the full 4' x 7' 
field of view (FOV) into the two subfields (channels 1 and 
2), each covering 4' x 3^5. We set the MOIRCS position 
angle to —90° so that the long side {— 7' side) of the 
FOV aligns to the east-west direction. In this setting, the 
channel 1 FOV is on the eastern (left) side (see Fig. lb). 

We used both the standard Xs-band filter and a narrow- 
band (NB) fiher called "CO" or "NB2288" . The NB filter 
has central wavelength Ac = 2.2881 /xm with a bandwidth 
of 0.023 /im. Note that these are the post-installation 
measured values (M. Akiyama, private communication). 
Individual NB exposures were 300 s. In total, we ob- 
tained 16 exposures, yielding a total integration time of 
4800 s. Shallow data were also acquired in the band 
to enable continuum subtraction of the NB data. The 
band exposures were taken in 55 s and 3-coadds. We ob- 
tained 9 frames in total (1485 s). We found that stars on 
three images were elongated due to poor guiding at high 
elevation (> 84 degrees). We rejected those frames during 
the final combine, resulting in a total integration of 990 s. 
We applied a standard circular-dither pattern with a di- 
ameter of 20". After the target observation, four position 
dithered data of the FS151 standard star were acquired in 
both the Kg and NB2288 bands. 

We note that 4C 23.56 is close (~ 1°) to a bright 
[K ^ mag) Mira-type star. Such bright stars are known 
to occasionally cause significant stray light on the im- 
age-'^. Fortunately, the effect was small during the night. 
However, stray light was seen in some of the data. It af- 
fected only a tiny portion (< 1%) of the total area, and 
we took care to mark the positions of any stray light be- 
fore identifying the emitters. Unfortunately, stray light 
caused a hole-like pattern near the radio galaxy in the Kg 
data. Due to this, the total magnitude for 4C 23.56 may 
be systematically wrong. For other targets, we confirmed 
that stray light does not affect the photometry. 

2.1.1. Wavelength Shift of NB filters 

In general, the central wavelength of NB filters will shift 
with the angle of incidence of the light to the filter surface. 
In our case, therefore, the observable wavelength window 
will change with the position of targets on the detector. 

In Figure 1 we show the simulated shift of the filter 
transmission curve with source position on the detector 
(the case for channel 2 is shown here: R. Suzuki, pri- 
vate comm.). We find that the redshift range of interest, 
2.47 < z < 2.49, is always covered by > 80% of the de- 
tector area. As a reference, the expected wavelength of 
the Ha emission from 4C 23.56 is plotted as the arrow in 
Fig.l. Note that the redshift of 4C 23.56 is 2.483 ±0.003 

^ See the Imaging Information page on MOIRCS website. 
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by Rottgering et al. (1997) from CiV, He II, and Cm] 
lines. The bottom of Fig.l shows the iso-redshift Imes 
of the observable Ha emission for the center of the filter 
transmission. The iso-redshift line for 2.483 goes across 
the center of the detector. In the upper panel of Fig.l 
we show that Ha emission at ±500 km offset from 
4C 23.56 can be observed everywhere but the very edge of 
the detector. The plot shows that the emitters blueshifted 
by ^ 1000 km s~^ relative to the radio galaxy will leave the 
filter at the center of the FOV, while for objects redshifted 
by ^ 1000 km s~^ escape our detection in the outer 1/4 
of the detector. Though MOIRCS uses two independent 
filters for the two channels, the difference of the center 
wavelength of these two NBs is only 3 A (~ 40 km s~^ 
difference at z = 2.49). It should not cause any systematics 
between the channels. 

We note that the velocity window of ±500 km s~^ is 
just adequate for cluster member selection. Simulations 
of X-ray clusters shows that the velocity dispersion of rich 
clusters generally gets smaller by about a factor of 2 as 
one goes from z = to 3 (e.g., Eke et al. 1998). Also, 
the previous Ha emitter survey for a dense protocluster at 
z = 2.16 by Kurk et al. (2004b) found a velocity dispersion 
of only 360 km s~^ (narrower than the width of the NB 
filter for their pre-selection), small for the inferred richness 
of the protocluster. Thus we would expect to detect a 
significant fraction of the protocluster Ha emitters within 
the MOIRCS narrowband window if its dynamical center 
is at the same redshift and located on 4C 23.56. 

2.2. MOIRCS Data Reduction by MCSRED 

Data reduction of the MOIRCS science frames was per- 
formed using the MCSRED package^, the IRAF"^ script 
pipeline for MOIRCS imaging data developed by the au- 
thor (IT). Here we briefly describe the procedures in the 
pipeline for the potential user's convenience. 

The first step is flatfielding the data. We use a self-flat 
(i.e., maky a flatfield image from the dithered object im- 
ages) for the Ks band. An object mask is essential for 
making good quality flatfields. To do this, we first make 
"quasi- flatfielded" images by subtracting the adjacent im- 
age in time from the image from which we make the object 
mask. Such a "quasi-flatfielded" image will show a flat sky 
and objects in the image as positive-negative pairs. Then 
the residual pattern in the sky of each quasi-flatfielded 
image is further subtracted by fitting a 3rd-order polyno- 
mial surface. This procedure is done independently for 
each quadrant, because MOIRCS is a 4-channcl readout 
system and any pattern related to the sky variation has 
discontinuities on quadrant boundaries. After replacing 
all the negative objects to the sky value, cosmic-ray re- 
moval is applied. Each resulting image is then smoothed 
by Gaussian kernel after a 3 by 3 boxcar median filter has 
been applied. At this stage, most objects are clearly seen 
on the smoothed image except very faint ones. Then any 
area that has a signal larger than the threshold of 1.25a 

^ Available at 

http: / / www.naoj .org/staff/ichi /MCSRED / mcsrcd.html 
^ http://iraf.noao.edu/ 



of the original pixel-to-pixel variation are replaced by the 
"mask" value of 1, with an additional expansion to the 
area (here we set 5 pixels) for large and bright objects. 
This quick object mask works well for the usual science 
exposures. We note that, though we can do the more accu- 
rate "second-pass" mask-making procedure after getting 
the final combined image, the result is almost the same 
and is not necessary for our data. 

After making the object masks, the flat frame is made 
from the science data itself with the mask applied. For the 
narrowband data we used dome flat data as we describe 
further below. After the flatfielding of the NB data, we 
further subtract the "median sky" image (made with the 
object masks) which is made from the subset (typically 
6 to 10) of the frames before and after the image for the 
sky subtraction. For the Kg data this procedure can be 
skipped because the sky pattern is relatively smooth (by 
the use of the self flat) and so we can subtract the sky 
in the next procedure reasonably well. Then we further 
subtract any residual sky pattern remaining on each im- 
age using a 3rd-order polynomial surface fitting for each 
quadrant. 

The resulting flatfielded, sky-subtracted images are 
then corrected to remove the geometric distortion. The 
distortion correction database, which is made from obser- 
vations of astrometric reference fields, is supplied by the 
Subaru Observatory. We then measure the position of 
some (typically 20 to 50) bright and compact objects on 
each image. The object catalogs for the images are fur- 
ther cross-matched with the catalog for the first image in 
the input list, and from the result we calculate the offset 
values we use to register the images to the first image. 
Finally, these images are combined, where we weight by 
the rms^ of the input images. 

A note on our NB2288 reduction: though we usually use 
the sky- or self-flats for the reduction, we used the dome 
flat data for correcting the pixel-to-pixel variation of the 
NB data. One reason is because the sky counts in the 
narrow-band filter are generally low compared with that 
of broad-band data. There is also the possibility that the 
sky may not be flat in the NB filter. Due to the shift of the 
observable wavelength window with position (see §2.1.1), 
the inclusion or exclusion of OH night-sky lines into the 
window becomes position-dependent. Although there is a 
small tilt pattern in the dome flats which is caused by the 
difference of the light path between dome data and object 
data, we estimate that the variation of the photometry 
due to the use of dome flat is < 3% (~ 0.033mag rms) by 
putting standard stars at several positions across the field 
of view. This is small compared to our required accuracy. 

2.3. Photometry and the Catalog 

The object detection and the photometry are done using 
the "double-image" mode of SExtractor 2.5.0 (Bertin & 
Arnouts 1996). For the detection image, we used the un- 
smoothed version of the narrowband data; therefore, our 
catalog is a simple narrowband-based selection. The Kg 
image is registered to this, and the seeing size is closely 
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matched (6.0 pixel, ^ 0.67 arcsec FWHM) for both im- 
ages. The object detection is done within the area that 
is covered by more than 12 frames (75%) for NB2288 as 
well as 4 frames (67%) for Kg. For the detection param- 
eters, we use the SExtractor DETECT_THRESH= 1.75 
and DETECT_MINAREA= 15 with a tophat detection 
filter (tophat_4.0_5x5.conv). 

The effective survey area is 23.6 arcmin^ which corre- 
sponds to a survey volume of 840.7 comoving Mpc^ at 
z=2.49 considering our NB2288 filter bandwidth. We used 
aperture magnitudes in a diameter of 12 pixels (= l."4 or 
11.5 kpc physical scale) for color measurement, and the 
MAG_AUTO values for derivation of the star-formation 
rates as well as for the stellar mass (§3.5). 

The standard star data for FS151 were reduced using 
the same flatfielding data we used for the object reduc- 
tions. The instrumental magnitude for FS151 is measured 
from the individual 4 frames without combining them be- 
cause the images are slightly defocused. These magnitudes 
were then averaged to get the final measurement. The in- 
strumental magnitude for each frame is measured by the 
"growth-curve" method which estimates the convergence 
magnitude with increasing aperture radius. We used the 
IRAF APPHOT package. 

The ifg-band magnitude for FS151 is taken from 
Leggett et al. (2006), with the -1-1.86 conversion from the 
Vega to AB-magnitude zeropoint applied (Tokunaga & 
Vacca 2005). For the NB2288 magnitude, we used the 
Kg magnitude and the low-dispersion template spectra of 
a G-3 dwarf star by WaUace & Hinkle (1997). As a fi- 
nal check, we compared the K-hand magnitude of some 
unsaturated stars in our catalog to the 2MASS photome- 
try. We find that the magnitudes are well matched at the 
■~ 0.01 magnitude level with a scatter of ~ 0.1. 

Good photometric error estimates are critical for ro- 
bust detection of emission-line objects at lower equivalent 
widths (EWs). In our analysis, we use the following two 
methods. First, we measure aperture photometry in 1000 
random (blank) regions to estimate the Icr level of the 
sky variation. This yields 5-a limiting magnitudes in a 
12-pixel diameter aperture of 22.3 and 22.9 for Ks and 
NB2288, respectively. Second, we estimate the scatter of 
the Xg— NB2288 color introduced by the photometric un- 
certainty. We generated 70000 fake objects with various 
magnitudes and placed them randomly in blank image re- 
gions. After rebinning the sample points in 0.5-magnitude 
intervals of NB magnitude, the scatter of the output color 
excess (Xg— NB2288) is evaluated by measuring the distri- 
bution's 68%(lcr), 95%(2cr), and 99.7%(3cr) levels directly 
from the simulation. The result of this simulation is used 
for the estimate of the error on the color excess and on 
the selection of emitters. 

2.4- Supplementary Data Used for Analysis 

In the current paper, we use some supplementary 
datascts (§3.1 & §3.4) and here we bricfiy describe them, 
though the full description will be presented elsewhere (I. 
Tanaka et al. in prep.). 



Deep J-band data were taken using the INGRID cam- 
era (Packham et al. 2003) on the William Herschel 4.2-m 
Telescope (WHT), La Palma. The data were acquired on 
the photometric nights of UT 2001 Sep 7th and 8th, un- 
der good seeing conditions (~ 0.6 arcsec). The INGRID 
camera has a HAWAII HgCdTe 1024 x 1024 array, cov- 
ering roughly 4' x 4' field of view (= 0." 2 38 /pixel) on the 
Cassegrain Focus of the WHT. The total integration time 
was 4800 s, resulting in a 50% detection completeness limit 
of 24.4 mag. The data were reduced in the standard man- 
ner, similar to the method for MOIRCS data reduction. 

Deep i3-band data were taken by the Suprime-cam on 
Subaru Telescope on the evening of UT 2009 Nov 17th 
under clear but a marginal seeing condition (0."9 ^ 1"2). 
It was taken as part of our Lya emitter search in this 
field. The net exposure time was 4500 s, having a limit- 
ing magnitude of 26.5 mag (5cr, 2."5-diameter aperture). 
The detail of the data reduction as well as the result of 
the broadband SED fitting analysis arc to be described 
elsewhere (I. Tanaka et al., in preparation). 

MOIRCS MOS spectroscopic observations of this field 
were done on UT 2006 September 7th (S06B-106: 
Kodama et al.); i.e., before the acquisition of the narrow- 
band data we describe here. Part of the target selection 
was done using a similar narrowband survey executed 
in September 2001 using the INGRID camera on WHT 
(IT and CP), which should sample the redshift range of 
2.429 < z < 2.489. The conditions were clear, with medium 
(~ 0"8) seeing conditions. We used the grism with a spec- 
tral resolution of ~ 1000 in the iiT-band window cover- 
ing the 2.0 — 2.4 fim spectral range. We used 0"8 wide 
slits with typical slit lengths of ^ 12". Each single expo- 
sure was 600 s to 1200 s, depending on the airglow line 
strength. The standard A-B dither pattern was applied 
during each exposure. In total 2-hrs of integration were 
gathered. The MOS data reduction was done in the stan- 
dard manner using the MCSMDP pipeline (pre-release 
version) by one of the MOIRCS builders (T. Yoshikawa, 
private comm.). Details of the MOS reduction proce- 
dure by MCSMDP can be found in the recent paper of 
Yoshikawa et al. (2010). 

Finally, the deep 24 /xm imaging from Spitzer/MlPS 
were obtained as part of GO program 30240 (PI: A. 
Stockton). The total exposure time in the central 2' x 2' 
of the mosaic is approximately 5 ks. The 171 separate ex- 
posures each had a background fitted and then subtracted 
before a median of the array was determined. This median 
was then subtracted from each exposure. The exposures 
were then mosaicked with MOPEX software"* with stan- 
dard inputs and resampling the pixels by a factor of 2. 

For the current analysis we focus on a subset of the 24 
/xm-selected sources that have clear IRAC counterparts. 
The MIPS photometry was extracted using DAOPHOT in 
IRAF to fit an empirical PSF model to multiple proximate 
(within 12") sources simultaneously. The MIPS source 
positions were pre-defined using the position of the closest 

http: / / ssc.spitzer.caltech.edu/dataanalysistools / tools / mopcx/ 
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IRAC counterpart. This photometry method mitigates 
the effects of source confusion. 

3. Results 

3.1. Emission-line Object Selection 

The selection of Ha emitters at z ~ 2.5 is based on ex- 
cess emission in NB2288 measured via i^s— NB2288 color. 
In Figure 2, we show the color-magnitude diagram for all 
objects detected in our observations. The overall distri- 
bution of the detected objects tends to become broader 
with increasing magnitudes due to the increasing photo- 
metric uncertainty. The intrinsic scatter is also seen for 
the brighter objects. The intrinsic scatter is caused by 
the variation of the spectral energy distributions (SEDs) 
in the Ks band. 

The two curves labeled as S = 2 & 3 in Fig. 2 are the 
conventionally used measures of the significance of nar- 
rowband excess (e.g., Geach et al. 2008). It is formulated 
as below, 

niK - mNB = -2.5log[l - S^iO-O-^l™^-™"^)]^ (1) 

where is the zeropoint of the narrowband image and S 
is the photometric uncertainty. The E = 2 line generally 
matches well with our simulated a = 2 curve. As the color 
scatter becomes large at NB2288> 22 mag, we limit the 
selection of emission-line objects to S > 2 and NB2288< 22 
(corresponding to 7cr level in 12-pixel aperture). We note 
that the use of our simulated a = 2 scattering curve docs 
not change the sample. At NB2288> 22, there are actually 
only a few emission- line object candidates with E > 2, due 
to the steep rise of the selection curve. 

At the bright end, the intrinsic scatter is relatively large 
due to the variation of the SEDs in the Kg band window 
and the selection of emitters with low-equivalent width 
may introduce spurious sources. We therefore introduce 
a further selection criterion using the rest-frame equiv- 
alent width, EWo > 50 A at z = 2.49. This corresponds 
to the Ks- NB2288 = 0.55 mag. The limiting Ha line 
flux of the emitter sample by these criteria is roughly 
7.5 X 10~^^erg s~^ cm~^. These values are computed us- 
ing the same formula used in Geach et al. (2008). 

We find 12 objects that satisfy the selection criteria 
set here. Removing one object which is found to be a 
diffraction spike by visual inspection, the remaining 11 
objects are cataloged as robust emission-line objects for 
further analysis. We note that the brightest object in the 
emitter sample is 4C 23.56. The observational properties 
of the emitter sample are summarized in Table 1. 

3.1.1. Spectroscopy 

Though spectroscopic observations for all the emitter 
candidates are not obtained yet, four of the 11 emission- 
line candidates have NIR spectroscopy. These data were 
taken before the acquisition of the current NB data. 
The one- and two-dimensional spectra for each target are 
shown in Figures 3 and 4. HAE#491 is the radio galaxy 
4C 23.56 itself. 

Due to the relatively low signal-to-noise ratio (SNR) of 
the spectroscopy, unambiguous identification of Ha via 



detection of the nearby [Nil] or [Sii] lines is only pos- 
sible for the radio galaxy. The objects HAE#526 and 
HAE#153 are relatively bright, and we can say this is not 
[OlIl]A5007 emission at z = 3.571 because of the lack of 
the [OlIl]A4959 near 2.267^m (as well as the lack of H/3 
at 2.222 /xm). Furthermore, the spectra of the latter ob- 
ject seems to have a hint of the broad component which 
would not appear if it is [O III] . The lines from these two 
sources may also be from the Paschen series. If it is the 
Paa (Pa/3) hnc, Pa/3 (Pa7) wiU be observed at 1.56 (1.95) 
fim. Unfortunately our spectra do not cover wavelengths 
below 2 fim. We have to make use of other information to 
exclude the latter possibility (see next section). For the 
object #354, the Ha signal at 2.288 iim is quite uncer- 
tain. However, we see a possible signal at 2.342 /im at the 
expected position of the [S II] line. Though both lines are 
quite low SNR, the inferred large [Sll]/(Ha-|-[Nll]) ratio 
of ~ 1 would indicate that the object harbours an AGN 
(Ostcrbrock 1989). We list the measured redshift values 
in Table 1 supposing the line is Ha. 

The velocity difference among these emitters is rela- 
tively small. They lie in the range of 2.287 — 2.289 iim 
(230 km s~^ range), which is in the red edge of the 
INGRID narrowband filter function wc used for pre- 
selection (2.25 — 2.29 /im). If these emission lines are not 
Ha but from the interlopers unrelated to 4C 23.56, we can 
expect a more scattered distribution in wavelength range 
across the INGRID filter function. This is a strong indi- 
cation that they are physically associated with each other 
and 4C 23.56. 

3.1.2. Broadband Colors of HAEs 

Our emission- line survey for Ha at z = 2.49 may con- 
tain interlopers at different redshifts. As we showed in 
the previous section, the spectroscopic detection of a sin- 
gle line remains ambiguous; it could be Paa for example. 
Geach et al. (2008) used color selection as well as pho- 
tometric redshifts to distinguish interlopers in their Ha 
emitter sample at z = 2.2. They estimated that, among 
their 55 "robust" sample galaxies, a third can be lines 
other than Ha (12 Paa, 6 Pa/3, and one [Olll]). Their 
result indicates that foreground objects (Paa objects at 
z = 0.22 and Pa/3 objects at z = 0.78) could be the major 
contaminants of our data. 

To make sure that our detected emission-line objects 
are mainly Ha emitters, we measure the continuum color 
of our emission-line candidates using our supplementary 
Subaru S-band and WHT J-band data (§2.4). The J- 
band data covers the central 3.'8 x 3.'4 area of the MOIRCS 
field of view (here we refer to the area as the "INGRID 
subfield" as shown in Figure 6), while the whole area is 
covered by B-band data. 

Though our selected emission-line candidates all lie in 
the INGRID subfield, 4 objects are not seen in our B- 
band data. Two more objects are overwhelmed by the 
PSF of bright nearby stars. We could measure a BJK 
color for five objects. Fortunately, all four spectroscopic 
targets are in this sub-sample. The B — J versus J — Ks 
two-color diagram is shown in Figure 5a. The Ha emitter 
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candidates are marked by filled boxes (red). 

As discussed in detail in the Appendix, the B — J versus 
J — Kg two-color diagram is found to be quite effective for 
rejecting objects at z < 1. Figure 5b is the BJK two-color 
distribution of galaxies in GOODS-North region from the 
MOIRCS Deep Survey (Kajisawa et al. 2009; Kajisawa et 
al. 2010, and references therein). The objects that satisfy 
(J — Ks) > 0.4 — J) — 0.2 are empirically shown to have 
spectroscopic redshift of z > 1 from Fig. 5b. 

The position of our 5 emission-line objects in Fig. 5a is 
quite favorable; with the exception of one object, they sat- 
isfy the BJK criterion. These 4 BJK emitters comprise 
our spectroscopic sample. Thus we can conclude that the 
emission line we detect spectroscopically is indeed Ha, 
not from the Paschen series. This is strong support for 
the existence of a protocluster at z = 2.5 in the 4C 23.56 
field. 

One object whose color is below the BJK line could 
be an interloper, though we cannot conclude this for cer- 
tain until its spectra is obtained. One such object can be 
expected from the field survey as we see in §4.1. During 
the following discussion we consider that all our selected 
emission-line objects are HAEs at z ^ 2.49, though a few 
of them might be foreground objects. 

3.2. Spatial Distribution of HAEs 

The spatial distribution of our HAE sample on the sky is 
displayed in Figure 6. The position of 4C 23.56 is marked 
with the large diamond. The distribution of our selected 
HAEs (shown with open red circles) is quite tightly clus- 
tered at the east to south-eastern side of the radio galaxy. 

The biased distribution of HAEs to the eastern (i.e., 
channel-1) side of 4C 23.56 is not due to the wavelength 
shift of the filter transmission window or any difference 
between channels. It is also not due to some unexpected 
trouble in the channel-2 filter, because 4C 23.56 and 
the adjacent two emitters are actually detected on the 
channel-2 side. The inhomogeneous distribution seems to 
be the refiection of the real distribution of the emitters. 

If we suppose that the majority of the protocluster 
members are strongly redshifted relative to 4C 23.56, the 
situation could change. As is discussed by §2.1.1, HAEs 
redshifted by >1000 km/s relative to 4C 23.56 will drop 
off from our detection window in the outer half of the 
FOV. In the inner half of the FOV we should detect ob- 
jects with large velocity difference. Though our three 
HAEs with spectroscopic redshift have small velocity off- 
set (< 230 km s^^), further spectroscopy for all the HAEs 
may uncover a biased velocity distribution. 

3.3. Offset Ha Emission in the "Eastern Clump" 

Here we comment on the objects in the area enclosed by 
the small white box in Fig. 6. There is a group of 4 emitters 
that lies within a 15"-diameter area (=300 kpc comoving), 
which is also shown in Figure 7 as a zoomed view. We call 
this area the "Eastern Clump" . The compact size of the 
HAE group suggests that they might eventually merge 
into one massive galaxy in the future and become (one 
of) the brightest cluster member(s). 



One notable character of HAEs in the Eastern Clump 
is that the light distribution of each HAE in the NB2288 
image is offset from the X-band object position, as shown 
in Fig. 7. Especially interesting is that the offset direction 
of emission-line gas in three HAEs on the west side of the 
clump in Fig. 7 is roughly aligned to the North to North 
West. 

To show these offsets more clearly, we ran SExtractor 
on each image independently and measured the positional 
offset of objects detected in both images. The result is 
shown in Figure 8. In the figure, the 1-, 2-, and 3-(t scatter, 
based on the objects that have similar magnitudes to the 
HAEs (dots), are shown by circles of dashed lines. The 
HAEs are marked as filled squares, with the objects in 
the Eastern Clump shown in Fig. 7 overlaid with asterisks. 
There are two more large-offset (> 3(t) objects seen in 
Fig. 8. One of them is indeed the object just south of the 
eastern clump. We confirm that the offset of the emission- 
line gas is not due to poor astrometric alignment of the 
two images; all the other objects adjacent to the Eastern 
Clump show the same scatter shown as dots here. 

The i^T-band light centroid can be considered the center- 
of-mass for each galaxy. Therefore, these offsets of the 
Ha emission indicate that the star-formation is occurring 
at the periphery of the galaxies. This could be due to 
interactions among the galaxies triggering star-formation 
(e.g., Sulentic et al. 2001; Cortese et al. 2006), though the 
apparent alignment of the offset direction is difficult to 
explain. The apparent offset of the Ha emission is roughly 
0"2 — 0"4, corresponding to a physical size of 1.6 — 3.2 kpc 
at z = 2.49. 

Another possible explanation may be that vast outflows 
of emission-line gas are being launched from these galax- 
ies. The direction of such outflows could be influenced by 
the direction that the group is moving towards the center 
of the cluster potential. If the cluster potential is already 
filled by dense intracluster gas, the expelled emission-line 
gas from these galaxies may leave a "wake" behind each 
galaxy. In the local universe, similarly offset light from 
emission-line gas is caused by ram-pressure stripping of 
galaxies by intracluster gas (e.g., Yoshida et al. 2002; Yagi 
et al. 2007; Sun, Donahue, & Voit 2007). For example, 
Cortese et al. (2007) discovered a trail of blue emission- 
line knots behind an object in Abcll2667 at z = 0.23. Our 
detected offset of the Ha emission in the Eastern Clump 
may be similar. 

A detection of extended X-ray emission from 4C 23.56 
intracluster gas is the key for falsifying the ram-pressure 
stripping scenario. The field was observed by XMM and 
reported by Johnson et al. (2007). Though they detected 
extended soft X-ray emission near the radio galaxy in a 
^ 30 ksec EPIC exposures, they concluded that inverse- 
Compton scattering of CMB photons is the most likely 
emission mechanism. Indeed the position and direction 
of the X-ray extension is well aligned with the radio lobe 
which extends ^ 0.5 Mpc north-east/south- west from the 
radio galaxy, and not with the direction where the Eastern 
Clump lies. For now, there is no evidence of thermal X-ray 
emission from hot intracluster gas. 
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3.4. HAEs and the Spitzer MIPS Data 

Spitzer MIPS 24 imaging has been used to probe 
star formation activity at z^2 (e.g., Daddi ct al. 2007). 
There is public, archived deep 24 /im MIPS data for the 
4C 23.56 field (§2.4). If the protocluster around 4C 23.56 
is characterized by active star-formation, we may also see 
it as an excess of faint MIPS sources. 

In Figure 9 we show our re-processed MIPS 24 /im im- 
age. The strong point source at the center of the image 
is 4C 23.56. We can clearly see an excess of faint MIPS 
sources around the radio galaxy compared to the outer 
parts of the image. We also find that 9 objects out of our 
11 HAE sample (shown as red open circles) have MIPS 
counterparts. However, the apparent excess of 24 /im 
sources does not seem to be due solely to the HAE coun- 
terparts. Here we examine the excess of MlPS-selected 
objects in the 4C 23.56 field. The 24 photometry and 
the inferred star-formation rates for the HAE counterparts 
are separately discussed in the next section (§3.5). 

To estimate the density excess of MIPS sources quanti- 
tatively, we have to consider the difference of the exposure 
coverage across the field. We first exclude any area that 
has exposure time less than 1/3 of the total integration. 
We then set the SExtractor object detection threshold to 
2.8(7. This is higher than usual so that the detection is 
not affected by the S/N difference due to the variation of 
exposure time across the field. To check the robustness 
of the detection threshold further, we reversed the pixel 
count of the MIPS image and executed the object detec- 
tion with the same SExtractor parameters. Besides a clear 
pixel defect, we detected no sources in the inverse image, 
indicating that the no fake object from the sky variation is 
included in the current analysis. Furthermore, we exclude 
all objects with flux density less than 90 /iJy, as they may 
suffer from the source confusion in crowded regions. 

We also introduce an upper flux density cut of 275 fiJy 
because most of our HAEs are fainter than 200 /iJy. Any 
object much brighter than 200 fj.3y has a high possibility 
of being in the foreground. The upper flux density cut is 
chosen to be roughly twice the median 24 fim flux density 
of our HAEs. 

Note that we apply the SExtractor detection directly 
to the 24 /im data for the source density excess estimate. 
For the SFR estimate from the MIPS photometry for our 
HAE sample we use DAOPHOT software as described in 
§2.4, because it provides a more robust flux estimate of 
each source. 

The resulting density map of the faint (90 — 275 ^Jy) 
MIPS sources is overlaid on Fig. 9 as contours. The distri- 
bution of sources has been smoothed by a Gaussian kernel 
of CT = 30". For reference, the position of the HAEs are 
shown as red open circles. It is quite clear that the MIPS 
source density excess closely matches the HAE distribu- 
tion. 

To quantify the density excess of MIPS sources further, 
we define the "excess region" as the 80"-radius region 
around the peak density contour, and compare the excess 
within the region with the "outer region" at larger radii. 



From the density of faint MIPS sources in the "outer re- 
gion" we estimate the expected number in the "excess 
region" as 15.3. On the other hand, we count 39 objects 
in the "excess region" , yielding a number enhancement of 
2.55lo;4g. This corresponds to a 4.5-ct excess by proper 
Poisson statistics (Gehrels 1986). As a crosscheck, we 
downloaded the deep MIPS data for GOODS-N/S fields^ 
and counted sources in the same manner, yielding the ex- 
pected count of 18.2 ± 0.7 for an area with the size of the 
"excess region". The source counts of deep MIPS data 
for ELAIS-Nl field by Chary et al. (2004) also shows the 
expected number of 17.0. Our number is consistent with 
those surveys, considering the uncertainty by the cosmic 
variation. 

The HAEs with MIPS counterparts (9) can only explain 
less than half of the total excess count 24) . This implies 
that our shallow Ha emitter survey may be missing an 
even larger number of the dust-enshrouded star-forming 
objects in the excess region. Future follow-up observations 
using ALMA etc would be quite interesting to understand 
these hidden dusty sources. 

3.5. Star- Formation Rate and Stellar Mass of HAEs 

In this section we estimate the star-formation rate 
(SFR) and the stellar mass of our HAE sample. We note, 
however, that the estimation of the SFRs and the stellar 
masses for our HAE sample are uncertain and should be 
treated with caution. 

3.5.1. Ha-based Star- Formation Rate of HAEs 

It is known that dust absorption lowers the measured 
SFR from Ha emission relative to the intrinsic SFR based 
on radio or FIR emission (e.g., Hopkins et al. 2001; Afonso 
et al. 2003; Erb et al. 2006; Daddi et al. 2007). The con- 
tamination of the [Nil] 6583A emission line to the narrow- 
band flux is also an unknown factor (here we set the factor 
as 0.2). The contribution of AGN to the ionizing radiation 
also adds uncertainty to the estimate. These uncertain- 
ties can be calibrated using future deep spectroscopy of 
complete HAE samples. 

Kennicutt (1998) gives the conversion formula from 
L{Ha) to the SFR as follows. 

SFR (Mq yr-i) = 7.9 x lO"''^ L{Ha) (erg s"^) (2) 

Though the conversion assumes the standard Salpeter 
(1955) Initial Mass Function (IMF), the IMF is known to 
overestimate the stellar mass. We here use the conversion 
correction for the Chabrier (2003) IMF by dividing the 
value by 1.7, for mainly practical reasons, so that our data 
can be compared to some recent Ha-based SFR estimates 
in the literature. Note that the dust extinction should 
affect the SFR estimate much more than the uncertainty 
due to [Nil] contamination. An AGN contribution can 
also exist for some galaxies; the unusually high derived 
SFR for 4C 23.56 is clearly due to the AGN contribution 
to the Ha flux. 

As for the extinction-correction factor, the BJK dia- 
gram (Fig. 5a) is a useful diagnostic for a part of our sam- 

^ available in http:/ /ssc. spitzer. caItech.edu/spitzermission/observingprograms/Iegac; 
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pie. In the figure we show the Ay = 1 extinction vector 
using the Calzetti (2001) extinction law. Assuming that 
the objects are pure disks with dust extinction, the color 
distribution of the HAEs indicates that the level of the ex- 
tinction ranges from to ^1. Using the Calzetti (2001) 
recipe and the conversion of Aj^q, = 0.82Ay , an extinction 
correction factor of 1 to 5 is estimated. 

For the remaining 6 objects we cannot estimate the ex- 
tinction. For comparison, the recent study by Onodera 
et al. (2010) estimated the level of extinction for their 
sBzA'-sclccted Ha emission-line objects that seem to have 
roughly similar stellar masses as our HAEs (§3.5.3). Their 
extinction estimate is based SED fitting with additional 
correction as prescribed by Cid Fernandes et al. (2005) 
and Savaglio et al. (2005). Their Table 4 indicates that the 
average of the extinction correction factor is 4.9. We ap- 
plied the correction factors of 4.9 for the objects we could 
not use the BJK diagram. The raw and the extinction- 
corrected SFR for each HAE is hsted in Table 2. 

3.5.2. Spitzer MIPS Star- Formation Rates 

The bolometric infrared luminosity of galaxies corre- 
lates very well with other measures of the star-formation 
rate and may be the most robust measurement of the SFR. 
The infrared light is emitted by dust-reprocessing of the 
emitted UV light from young stars and therefore peaks at 
rest-frame wavelengths between 100 and 300 fim. In lieu of 
observations at these frequencies, the longest wavelength 
datapoints need to be extrapolated using assumed SED 
models. For this radio galaxy field we have used the ob- 
served 24 /im imaging to infer the SFRs for the discovered 
Ha emitters. 

We used the SED template for Arp 220, a possible 
low-redshift analogue, to correct the 24 fim to rest-frame 
8.0 /im luminosities (Rieke et al. 2009). From these val- 
ues we estimated the bolometric infrared luminosity using 
the Chary & Elbaz (2001) and Dale & Helou (2002) SED 
models, following the prescription by Daddi et al. (2007)^. 
The resulted infrared luminosities from the two SED mod- 
els are averaged. Finally, we use the Kennicutt (1998) Lir 
to SFR conversion below, 

SFR (Mo yr-i) = 3.0 x lO'^"^ ipiR (erg s"^). (3) 

Note that the original formula is divided by 1.5 for 
Chabrier (2003) IMF (Yoshikawa et al. 2010). These con- 
versions and extrapolation result in SFRs which are uncer- 
tain by up to factors of a few (Papovich et al. 2007). For 
our purposes, to estimate the FIR output from likely pro- 
tocluster galaxies, these uncertainties are tolerable. The 
Herschel Space Telescope has already quantified these un- 
certainties for other fields (e.g., Rodighiero et al. 2010) 
and may eventually provide longer wavelength data for 
these and other protocluster galaxies. The result of the 
MlPS-based SFR is also listed in Table 2. 

The resulting MlPS-based SFRs are significantly (a fac- 
tor 2-14 with median of 3.6) higher than the SFRs based 

^ The formula (4) in Daddi ct al. (2007) has a possible typo. The 
constant should be -1-1.745, not -1-1.23, to match the desciption 
in the text and to reproduce the relation in their Fig. 7 



on extinciton-corrected Ha. We will discuss this further 
in §4.3. 

3.5.3. Stellar Mass of HAEs 

Next we estimate the stellar mass of our HAE emit- 
ters. We employ the conversion formula for the stellar 
mass from the Xg-band photometry by Daddi ct al. (2004). 
They showed that for galaxies at 1.4 < z < 2.5 the stellar 
mass is reasonably estimated using the formula below, 

log(M,/10" Mq) =-0.4 (K'^%,~Kn) 

+0.218 [(z-/Oab -2.29] (4) 

where Kn is the X-band magnitude (Vega scale) of the 
average 10^^ M0-mass objects of their sample (=20.14 
for the average redshift of (z) = 1.9). Considering the 
difference of the luminosity distance between (z) = 1.9 and 
2.49, and the average distribution of {z — A')ab color of 
1 from the Fig.7 of Daddi et al. (2004), we derive the 
estimate of the stellar mass at z = 2.49 as 

log(M,/10" M0)2.49 = - 0.4 + 8.8 (5) 

where is in AB-magnitude scale. 

As a check of the validity of the formula, we com- 
pare it to the stellar mass estimates based on SED fitting 
the 12-band photometry for the GOODS-North MOIRCS 
Deep Survey Catalog (for details of the model fitting, see 
Kajisawa et al. 2009; Kajisawa et al. 2010). Figure 10 
shows the stellar mass versus if -band magnitude of the 
cataloged objects in the redshift range 2.3 < z < 2.6, for 
objects with photometric redshifts (open boxes) and spec- 
troscopic redshifts (crosses). Though the Daddi estimator 
by equation (4) fits the overall distribution of the photo- 
metric redshift sample reasonably well, the slope seems 
a bit shallower than the data. Therefore we fit all the 
photo- z-bascd data in Fig. 10 and got the best fit line as 

log(M,/10" Mq)2.49 = -0.5127 Al°]^-h 11.312. (6) 

We use this formula for the conversion from A'-band mag- 
nitude to total stellar mass. The rms scatter of the fit is 
calculated and found to be only 0.2 dex, consistent with 
the value quoted by Daddi et al. (2004) (but note that the 
scatter increases to 2.4 dex if we use brighter subsaniple of 
Ks < 23). The resulting stellar mass for each HAE galaxy 
is listed in Table 2. We note that the ATg-band magni- 
tudes for HAEs will be affected by the Ha emission. The 
rest-frame equivalent width of 100 A will affect the mass 
estimate by 11%. The median value of the offset due to 
line emission is ^ 20%. This is small compared to the 
mass estimation uncertainty of > 0.2 dex. We also calcu- 
late the specific star-formation rate (SSFR) for the HAEs. 
The results are shown in Table 2. 

4. Discussion 

4.1. Comparison to Field Counts 

Our discovery of 11 HAE candidates to a limiting fiux 
of 7.5 X 10~^^ erg s^^ cm~^ in 23.6 arcmin^ area should be 
compared with the number expected in the general field. 
Geach et al. (2008) provide us with excellent reference 
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data. Their HiZELS survey measured the number density 
of Ha emitters at z = 2.23 over 0.6 deg^ using the narrow- 
band filter (H2SI, Ac = 2.121 ^im, SX = 0.021 ^im) to a 
line flux limit of 10~^^ erg cm~^. Their survey field 
is 91.5 times wider than our own and their filter has a 
similar bandwidth to our NB2288 filter. They detected 
180 emitters including non-Ha (e.g., Paschen series, [Fell], 
[Olll]) objects under their selection criteria of observed 
EW of 50 A (= 14.3 A in rest frame at z = 2.5). As our 
data is deeper than theirs, we set our detection limit to 
match theirs. To a NB magnitude of 21.3 mag we find 5 
emitters. Next we have to estimate the how many of their 
180 emitters would satisfy our threshold of rest-frame EW 
of 50 A. In the filter bandwidth they used, the rest-frame 
EWo of 50 A (at z = 2.49) corresponds to a narrowband 
excess of 0.56. We estimate the number of emitters that 
satisfy the EWq > 50 A to be about 90 from their Figure 1. 
Scaling the number to match our field of view, we derive 
the expected number of the field emitters in our field to be 
~ 1. Note that, similar to our sample of HAE candidates, 
this field count includes interlopers too. Considering our 
total of 5 emitters, we can say that our field contains 5 
times the expected field count. 

Recently Geach et al. (2010) have published more gen- 
eralized HAE counts from their HiZELS and other sur- 
veys. Their Table 2 gives the expected field HAE number 
counts at the limit of 1 lO-^*^ erg s~i cm"^ (^-^Wo > 10 A) 
to z = 2.2. Extrapolating their result to z = 2.5 by third- 
order polynominal and applying our survey size, we get 
the exptected field count in our survey area as 1.2to;46- 
If we again suppose roughly half of them can be observed 
as HAE in our EW cut of EWo > 50 A, the expected field 
count could be only 0.6. Thus the result suggests the 
overdensity even more strongly. 

A deep and wide narrowband HAE survey using the 
new instrument HAWK-I was recently published (Hayes, 
Schacrer, & Ostlin 2010). A field of 7.'5 x 7.'5 in GOODS 
South was surveyed using the narrowband filter NB2090 
(Ac = 2.095 /im; SX = 0.019 ^m). Their data are extremely 
deep, with a 5ct limiting magnitude of 24.6 AB correspond- 
ing to a line flux limit of 6.8 x 10""'^^ erg s~^cm~^, nearly 
an order of magnitude deeper than our own data. The 
inset histogram in their Figure 1 shows that the number 
of emitters they detected at NB< 22 is 9, under their EW 
cut of EWo > 20 A (corresponding to the emission-line 
excess of 0.29). In their filter bandwidth our selection 
criterion of EWo > 50 A corresponds to an emission-line 
excess of 0.60 mag. We estimate that 5 of their emitters 
satisfy our selection criterion of EWo > 50 A from their 
Figure 1. Considering the difference in survey areas, the 
expected field count of the NB emitters for our data is 
estimated to be 2.5 ± 1. This estimate is consistent with 
the results from another deep survey, carried out on a 
smaller field, by Moorwood et al. (2000). Our count of 
11 HAEs is therefore 5.3^2'g times (or 3ct) more than the 
field count, which is based on the proper treatment of the 
small number Poisson statistics by Gehrels (1986). 

We note that the redshift evolution of the HAE number 
count suggested by Geach et al. (2010) is quite strong: 



extrapolating from redshifts below z = 2.2 to z = 2.5, the 
predicted number counts of HAEs with flux larger than 
1 10~^^ erg s~^ cm~^ decreases by 60%, from 8582 to 
^5380 deg~^z~^. Although the evolution in number den- 
sity depends on the limiting fiux, this correction for the 
difference in redshifts (from z = 2.2 to 2.5) will also reduce 
our field count estimate based on Hayes et al. (2010). Our 
estimate of the overdensity of a factor 5 should therefore 
be considered conservative. 

As we show in Table 2, about half of our HAEs seem to 
have a stellar mass of ^ 10^^ M0. Such massive objects at 
z ~ 2.5 are relatively rare in the field. For example, if we 
count such objects using the MODS catalog (§3.5.3), we 
find only 10 in their 103-arcmin^ survey area at 2.3 < z < 
2.6 (Fig. 10). By scaling the value to our survey area of 
23.6 arcmin^ and 2.47 < z < 2.49, we estimate the expected 
number in our FOV to be 0.15. The excess of such massive 
galaxies in the 4C 23.56 field is clear. Recalling that these 
massive galaxies in the 4C 23.56 field are Ha-selected, 
there could be more inactive massive galaxies missed by 
our search. We are planning to execute multi-wavelength 
observations of the 4C 23.56 field. The data will help 
identify such inactive massive galaxies, and will reveal the 
whole picture of the protocluster. 

4-.2. Comparison with Other Protoclusters 

Although there are several known protoclusters at z > 2, 
there is only one other field that has been searched for 
HAEs. Kurk et al. (2004a) searched for HAEs in their 
discovered protocluster surrounding PKS 1138-262 at z = 
2.16. Here we compare the number of detected HAEs. 

According to their Figure 6 and Table 5, their selec- 
tion of EWo > 50 A and S > 3 combined with a mag- 
nitude limit of NBvcga < 20.2 would match the selection 
we use here. There are ^ 10 HAEs meeting these criteria 
in the PKS 1138-262 field, roughly the same number as 
the 4C 23.56 field. Their HAE survey field is roughly haff 
of ours (but wider than our "INGRID subfield" where all 
our HAEs lie). However, if we use the distribution of their 
Lya emitters as a proxy for the HAE distribution (their 
Fig. 10) there would be at most an additional two HAEs. 
This implies that the richnesses of the 4C 23.56 and the 
PKS 1138-262 fields are likely comparable. 

There are differences in the properties of the HAEs in 
the PKS 1138-262 and the 4C 23.56 protoclusters, pri- 
marily the spatial distribution of the HAEs relative to 
the radio galaxies. In the PKS 1138-262 system the 
HAEs are rather homogeneously distributed around the 
radio galaxy, with some evidence for higher surface den- 
sity closer to it. In the 4C 23.56 field the radio galaxy lies 
at the western edge of the HAE distribution. This does 
not appear to be a result of observational bias as described 
in §2.1.1 and §3.2. 

This may imply a different evolutionary future of 
the host galaxy of each radio source. For exam- 
ple, the PKS 1138-262 host may evolve into the 
centrally-dominant cD-like galaxy relatively early while 
the 4C 23.56 host may remain one of several dominant 
galaxies until they merge into a single larger galaxy at 
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the center of the potential The host of 4C 23.56 is 
indeed ~ 1.5 magnitudes less luminous in Kg than the 
PKS 1138 host {Kg ~ 18.66 after removing nuclear com- 
ponent and emission-line contribution; Pentericci et al. 
1997). Furthermore, the Kg = 18.66 magnitude for the 
host of PKS 1138 is unusual; its apparent magnitude is as 
bright as the brightest cluster galaxies at z = 1.5 (Collins et 
al. 2009). If we use formula (6) in §3.5.3, the stellar mass 
of the 4C 23.56 and the PKS 1138 hosts are estimated to 
be 4.7 X IO^^Mq and 5.5 x IO^^Mq, respectively. For com- 
parison, Seymour et al. (2007) estimated the stellar mass 
of both host galaxies based on their Spitzer multi-band 
photometry, finding 3.9 x lO^Mg and < 1.3 x 10^^^^^ 
for 4C 23.56 and PKS 1138 respectively. 

Another difference is in the number of HAEs with strong 
(EWo > 100 A) emission. In the 4C 23.56 field, 10 of 11 
HAEs have rest-frame EW^ > 100 A. In the PKS 1138 
field the number of objects that satisfy EWo > 100 A and 
NBvoga< 20.2 is 5. As the equivalent width relates to the 
Ha luminosity, the excess number of strong EW objects 
in the 4C 23.56 field implies more active star formation 
occurring in that protocluster. We will address this quite 
interesting topic more in the next section. 

Recently, Hayashi et al. (2010) reported an excess of 
[O II] emitters in the field of the distant X-ray cluster 
XMMXCS J2215. 9-1738 at z = 1.46. The distribution of 
[O II] emitters are prevalent from the outskirt to the core 
region of the X-ray cluster. They showed that the frac- 
tion of the [O II] emitters among all cluster members as 
well as the level of the star-formation are shown to be 
both nearly constant as a function of clustercentric radius 
(see also, Tran et al. 2010). 

Their Figs. 5-7 indicate that the distribution of the 
emitters well traces the overall distribution of all clus- 
ter members. The distribution of Ha emitters in Kurk 
et al. (2004a) also closely matches the photo-z-selected 
galaxy distribution as shown by their Fig. 6 (Tanaka et 
al. 2010). We note, however, that the case of the cluster 
XMMXCS J2215.9-1738 may not be the representative as 
a cluster at z 1.5. Hilton et al. (2010) has shown that 
the system has a possible bimodal kinematic structure, 
with the position of the cluster BCG significantly offseted 
from the center of the thermal hot X-ray gas. And also, 
the NB HAE survey for another cluster XMMU J2235.3- 
2557 at z 1.4 by Bauer et al. (2010) has show that the 
HAE distribution just avoids the central ~ 200 kpc of the 
cluster core (see also Lidman et al. 2008; Strazzullo et al. 
2010). Although we are still not sure about whether the 
HAE can trace the overall distribution of cluster galax- 
ies from the core to the outskirts at 2 > 2 universe, both 
studies could consider the supporting case that the distri- 
bution of HAEs in the 4C 23.56 protocluster also traces 
the underlying galaxy distribution. Future analysis based 
on the broadband imaging will enable us to address the 
issue. 

Hayashi et al. (2010) contrast the activity of galaxies 
in the cluster core region by comparing to similar work 
for a z = 0.8 cluster by Koyama et al. (2010). At even 



lower redshift, Kodama et al. (2004) studied the HAE 
distribution for CI 0024-1-1652 at z = 0.4, showing quite 
strong suppression of activity in the cluster core region. 
It is quite interesting to connect the results by these works 
from z = 0.4 (Kodama et al. 2004), z = 0.8 (Koyama et al. 
2010), z= 1.4- 1.5 (Hayashi et al. 2010; Bauer et al. 2010), 
z = 2.1 - 2.2 (Kurk et al. 2004a; Gobat et al. 2010), to 
our z = 2.5 studies as a snapshot of the cluster evolution 
in each epoch (Tadaki et al. 2010). The implied rise of 
star-formation activity in the cluster core with redshift 
and its variety suggests that, while it may indicate that 
we are gradually approaching the era of galaxy formation 
for members in the cluster core (Kodama & Bower 2001), 
it is important to evaluate the physical state (richness, 
dynamics, mass etc.) for each protoclusters in detail to 
draw a general pictures of early history of galaxy evolution 
in cluster core. The detailed studies of larger samples are 
clearly needed to corroborate this apparent evolutionary 
trend. 

4.3. Star- Formation Activity of HAEs 

In §3.5 we calculate the SFR and the stellar mass of 
our HAE sample. Although we admit that the estimate of 
these values is rather crude, here we try to characterize the 
star-formation of the HAEs in 4C 23.56 field by comparing 
them to the other studies. As is discussed in §3.5, the 
main uncertainty for Ha-based SFR is in the estimate of 
the dust extinction, as indicated by the MIPS detection. 
The unknown fraction of the AGN contribution to our 
HAE sample (it can be non-negligible: see Lehmer et al. 
2009; Digby-North et al. 2010; Lemaux et al. 2010) is also 
an issue we need to address by the future observation. 

In Figure 11a, we plot the estimated star- formation 
rate of our HAEs against their stellar mass. The mini- 
mum SFR in our sample is 15 M0 yr^^, with a median 
value of 45 (182 for extinction corrected, 746 for MIPS) 
Mq yT~^. In the figure, the brightest object with the 
highest SFR, is 4C 23.56 (note that the mass estimate for 
this object is probably overestimated, not only due to the 
inclusion of AGN but also due to the effect of the poor 
sky subtraction: see §2.1). The line on the figure is the 
average mass-SFR relation for z '--^ 2 BzK-sclcctcd galaxies 
with MIPS detections by Daddi et al. (2007). The SFRs 
based on both indicators show larger values than expected 
from the z ~ 2 mass-SFR relation by Daddi et al. (2007). 
Especially, the MlPS-based SFR shows large excess. This 
trend seems real, because our MlPS-based SFR is derived 
from the same equation as used by Daddi et al. (2007). 
This may indicate that some galaxies in the protocluster 
4C 23.56 have higher star-formation activity than the field 
counterparts with similar mass. 

The observed large excess of the 24^m-based SFR rela- 
tive to Ha-hased one may indicate that the star-formation 
activity in our HAE sample is largely (or selectively) hid- 
den by dust. At lower redshifts such hidden star- formation 
activity revealed by infrared observation is often reported 
(e.g., Valdes et al. 2005; Geach et al. 2006; Garcia-Marfn 
et al. 2009; Garn et al. 2010). On the other hand, the 
recent calibration of the MIPS 24/xm-based SFR for star- 
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forming galaxies at z ~ 2 by the Herschel PACS observa- 
tion has revealed that the SFRs from 24/zm alone could 
be overestimated by a factor of ~ 2 — 7.5 (Nordon et al. 
2010; Rodighicro et al. 2010). Our observed trend can be 
due to the large systematics in the estimate of the total 
infrared luminosity from the rest-frame 8fim photometry. 
For now we do not know whether the HAEs in the 4C 23.56 
field harbours the star-formation significantly hidden by 
dust. The direct far- infrared observation by ALMA will 
address the issue. 

In Fig. lib wc show the L(HQ;)-based specific star- 
formation rate (SSFR) as a function of stellar mass. For 
comparison, we plot the SSFR values for field Ha emit- 
ters from Forster-Schreiber et al. (2009) and Onodera et 
al. (2010). The Forster-Schreiber et al. (2009) points 
include two types of the extinction correction, which is 
shown as connected points with vertical lines in the figure, 
while Onodera et al. (2010) used the correction (described 
above) which could be considered the middle between the 
two extreme cases from Forster-Schreiber et al. (2009). 

The measured relation between SSFR and mass devi- 
ates clearly from that observed in similar studies of low 
redshift clusters. For example, Patel et al. (2009) mea- 
sured the SFR and SSFR of a rich cluster at z '--^ 0.8 from 
MlPS-based photometry, and confirmed that the SSFR of 
massive (log M/Mq > 10.8) galaxies in high density re- 
gions is as low as for field objects with similar mass in 
local universe. Also at z > 1, the SSFR estimate of the 
member galaxies in the irregular cluster CL 0332-2742 at 
z ~ 1.6 by Kurk et al. (2009) show that significant member 
galaxies have much lower SSFR than field counterparts at 
similar redshifts, despite that the system is quite irregu- 
lar and far from a virializcd structure. Recently, Tanaka 
et al. (2010) studied galaxies in the protocluster around 
PKS 1138 at z = 2.16 and showed a similar trend of the 
suppressed star-formation activity in the cluster compared 
to the field environment (see also, Zirm et al. 2008). These 
studies may indicate that the star-formation activity in 
cluster environment is suppressed up to z '-.^ 2. 

To avoid the uncertainties introduced by estimating 
the SFR and the stellar mass, we also compare the star- 
formation rates directly using the observed values. In 
Figure 12 we show the observed if-band magnitude versus 
the Ha luminosity (L(Ha)) of the same dataset as shown 
in Fig. lib, but limited to objects at z > 2. Here the data 
from Kurk et al. (2004a) is also added for comparison. 
Within the uncertainties, all HAEs in the 4C 23.56 field 
are clearly as luminous as the brightest Ha sample in the 
general field results, again suggesting an excess of active 
star formation in the 4C 23.56 field. 

Perhaps more interestingly, in Fig. 12 we compare the 
HAEs in the 4C 23.56 field with the PKS 1138 field data. 
The majority of HAEs in the PKS 1138 protocluster have 
lower Ha luminosities than both field objects and the 
HAEs in 4C 23.56 despite being within a similar K mag- 
nitude range. It indicates that the HAEs in the PKS 1138 
field show suppressed star formation. This is actually the 
same conclusion reached by Tanaka et al. (2010), but their 
work is based on broadband SED fitting, i.e., from an in- 



dependent measure. 

The cosmic time between the two protoclusters is 
only 0.4 Gyr. This is actually long enough to dis- 
cern changes of the star-formation activity if their star- 
formation timescale is short. Tanaka et al. (2010) sug- 
gest that the star-formation timescale for galaxies in the 
PKS 1138 field is much shorter than for galaxies with sim- 
ilar mass in the GOODS fields. Their Figure 13 demon- 
strates the rapid decline of the SFR in the protocluster 
compared to that of GOODS, and shows that the time 
difference of 0.4 Gyr is long enough to observe the de- 
cline. For example, if we suppose that these star-forming 
galaxies in protoclusters are formed at z 3, the increase 
of the star-formation rate between z = 2.16 to 2.5 (0.4 
Gyr) is close to an order of magnitude. 

This argument could of course be too simplified, as it as- 
sumes that any galaxies in the two clusters have almost the 
same formation period around z 3. The actual forma- 
tion of massive galaxies in the cluster environment could 
be much more complicated, and it could be coupled with 
the fundamental cluster properties (mass, merger history, 
dynamical state, and so on) . As we already discussed, the 
difference of the star-formation activity between the two 
protoclusters might be related to the different inherent 
cluster characteristics. The PKS 1138 host is exception- 
ally massive and already lies near the apparent center of 
the galaxy distribution (Kurk et al. 2004a; Seymour et 
al. 2007), while the 4C 23.56 host lies at the edge of the 
HAE (and presumably overall galaxy) distribution (§3.2 
& §4.2). The soft X-ray emission from the hot intraclus- 
ter gas is (marginally) detected in PKS 1138 (CariUi et al. 
2002), while for 4C 23.56 it is not (Johnson et al. 2007). 
The future of the PKS 1138 protocluster may well be an 
X-ray bright, well evolved rich cluster with a massive cD 
galaxy lying in the center, while that of 4C 23.56 might 
be an X-ray faint cluster with a few dominant giant ellip- 
ticals near the center of the potential, such as the Coma 
or Virgo clusters. 

5. Summary 

We have discovered a significant excess of near- 
infrared selected emission-line objects around the ra- 
dio galaxy 4C 23.56 at z = 2.486. We detect 11 ro- 
bust emission-line galaxies with a line flux greater than 
7.5 X lO^^^crg s^^ cm^^ in our MOIRCS imaging obser- 
vation of the 4' X 7' held around the radio galaxy. We 
show that this is a significantly higher (5 times more) de- 
tection rate compared to other near-infrared emission-line 
surveys in non-radio-galaxy fields. Three of the emission- 
line galaxies are found to have a spectroscopic redshift of 
z = 2.49. We therefore conclude there is a protocluster 
physically associated with 4C 23.56. 

The distribution of HAEs on the sky is confined to a 
relatively small area 1.2Mpc-radius region) near the ra- 
dio galaxy, with 4C 23.56 placed at the western edge. We 
show that this is not due to observational bias, and con- 
clude that the protocluster is physically associated with 
the radio galaxy. The offset position of the radio galaxy 
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is interesting, because the hosts of powerful radio galax- 
ies are usually very massive and often assumed to be the 
dominant object (e.g., cD galaxy) in the center of the 
cluster. We also find that there is an interesting compact 
group of four HAEs, all the members of which show offset 
Ha emission relative to the isTs-band light position. This 
may indicate strong star-formation on the edge of each 
galaxy which is occurring due to strong interactions with 
the group. A more exotic interpretation might be that 
there is ram-pressure stripping of the galactic-scale Ha- 
emitting gas by the hot intraclustcr gas. Further follow-up 
is necessary to understand the phenomenon. 

Our analysis of deep Spitzer MIPS 24 yum data for 
4C 23.56 shows that all but two of the HAEs have 24 //m 
counterparts. The inferred star-formation rate (SFR) 
based on the MIPS 24 jum photometry is quite high (me- 
dian of > 700 M0 yr~^), indicating that intense starburst 
activity is going on in each HAE. The MlPS-based SFR is 
on average 3.5 times higher than the value inferred from 
the observed Ha luminosity. The strong effect of "selec- 
tive" dust extinction is evident for these galaxies, imply- 
ing the metal-rich and evolved populations for HAEs, as 
suggested by Kurk et al. (2004a). 

The analysis of the spatial distribution of the MIPS 
sources has also revealed that there exists a strong excess 
around 4C 23.56. Notably, the distribution of sources that 
have 24 fim fluxes as faint as the HAEs (90 — 275 /iJy) 
shows the same offset distribution relative to the radio 
galaxy, and the position of the density peak nicely coin- 
cides with the HAE distribution. This is another indica- 
tion that the strong star-formation activity is going on in 
the protocluster. 

The comparison of the star-formation rates of HAEs in 
the protocluster 4C 23.56 at z = 2.48 to those in another 
protocluster PKS 1138-262 at z = 2.16 and to those in 
the general field at z ^ 2 shows that: (1) HAEs in the 
4C 23.56 field are as active as those in the general field, 
and that (2) HAEs in the protocluster PKS 1138-262 show 
a clear deficit of objects with high Ha luminosities. The 
clear difference in star forming activities between the two 
protoclustcrs at similar rcdshifts may indicate that we are 
stepping into the epoch of major star formation in massive 
cluster galaxies. Given the fact that there is a small but 
significant time difference of 0.4 Gyr between the two red- 
shifts, such variation in star formation activity may be in- 
terpreted as a time evolution where star formation activi- 
ties to form massive galaxies can decay rather sharply dur- 
ing this epoch (2 < z < 3) in proto-clustcrs (e.g., Kodama 
et al. 2007; Tanaka et al. 2010). Alternatively, we may be 
just seeing the intrinsic scatter of cluster properties at this 
high redshift where the onset of cluster formation and the 
subsequent evolution depend on the mass of the systems 
and/or structures in the surrounding environment. Such 
dependence is actually inferred in lower- redshift clusters 
at z < 1 (e.g., Urquhart et al. 2010). However, we do not 
know yet the global properties of these protoclustcrs, such 
as total number of member galaxies, cluster masses, and 
the dynamical states, nor do we yet have a statistical sam- 
ple of protoclustcrs at this interesting redshift regime to 



disentangle these two effects. We are desperate to obtain 
such information by extensive near-infrared spectroscopy 
and to increase the sample of proto-clusters at z > 2. 
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Appendix 1. The BJK Color Selection for 
Rejecting Objects at z < 1 

In order to reject the possibility that our single line 
detection by NB filter or the spectra is mainly from the 
low-redshift Paschen series interlopers, we introduce the 
broadband color selection by the B — J versus J — Kg two- 
color diagram. Here we show its effectiveness using the 
large photometric/spectroscopic catalog for the GOODS- 
North region. 

We make use of the MOIRCS Deep Survey Catalog 
which contains 13-band photometric data as well as 
the extensive compilation of spectroscopic redshift data 
(Kajisawa private com.; Kajisawa et al. 2010). In 
Figure 5a, we plot the color of all galaxies with Kg < 23 in 
GOODS-N (gray crosses). The galaxies at 0.2 < z < 0.25 
and 0.75 < z < 0.8 (spectroscopically measured), each of 
which includes the Paa and Pa/3 interlopers, are shown 
by blue filled triangles, while the objects at 2.3 < z < 2.6, 
which include the 4C 23.56 redshift, are shown by red 
filled circles. There is a clear line between the two sam- 
ples. Here we define it as, 

(J-Xs) = 0.4 (B- J)-0.2, (Al) 

which is shown as the magenta dotted line in the fig- 
ure. The insetted diagram in Fig. 5b shows the fraction of 
galaxies whose J — Kg color is redder than the line defied 
here {'^BJK" galaxies) as a function of redshifts (spectro- 
scopic redshift sample only). Almost no objects at z < 1.2 
can be BJK galaxies, while more than 80% of objects at 
z ~ 2.5 satisfy the criterion. 

We also calculate the behavior of galaxy colors in the 
BJK plane using the population synthesis SED models 
by Kodama & Arimoto (1997) ovcrplotted in Fig. 5a. In 
the figure, each model track indicates, from bluer to red- 
der, the "pure disk" model with continuous SFR (bluest), 
"50% disk -I- 50% bulge", "20% disk -t- 80% bulge", and 
"pure bulge" model with the formation redshift of 5 (red- 
dest, and mostly out of the figure). In each model, sym- 
bols are put from z = 0.1 to 4, with the interval of 0.5 
from z = 0.5 (each line starts from z = 0.1 with a symbol, 
and ends at 4.5 without symbol). A part of each model 
track has the region with cyan-colored solid line, which 
indicates the region at 2 < z < 3.8. The extinction vector 
for Ay = 1.0 based on Calzetti (2001) is also plotted. 

The behavior of the model tracks also confirms that the 
our proposed dividing line for BJK galaxies works reject- 
ing objects at z < 1. The model tracks at 2.3 < z < 2.6 
(the same redshift range as red filled circles in Fig. 5b) 
tend to go a bit further from the BJK line than the ac- 
tual distribution of objects at the same redshift range in 
Fig. 5b. This may be explained, though, by a difference in 
the metallicity between the modeled and observed galax- 
ies. More advanced modeling is beyond the scope of this 
paper, however. 



Appendix 2. Objects Identified in Previous 
Studies 

As is described in §1, the 4C 23.56 field was also ob- 
served previously by KC97, Kajisawa et al. (2006b), and 
Stockton et al. (2004). Here we briefiy comments on the 
objects detected in each study. 

A. 2.1. Objects in KC97 

KC97 studied the 1.25 arcmin^ region around the radio 
galaxy using deep U'VIK imaging data from 2.2- & 3.6-m 
telescopes. They also searched for the Lya emitters using 
the narrowband filter at 4229A on the UH 2.2-m telescope. 
They cataloged 9 possible Lya-emitting companions and 
15 EROs selected hyl — Ks colors around the radio galaxy. 

We first looked for the counterparts of their Lya emit- 
ters in our K^ image. However, we can only identify one 
object (ID#110 in their catalog). The K^-NB color of 
the object is 0.14, so no clear sign of NB excess is seen in 
Figure 13 (red filled box). 

Next we identify their EROs. Of their 15 objects, 6 ob- 
jects are not seen at all in our I^s nor NB images, and 5 ob- 
jects are quite faint and subsequently not in our SExtrctor 
catalog. Note that their data is extremely deep, with a 3(j 
limiting magnitude oi K = 24.3 mag, while our data only 
reaches 23.5 mag. Among the 4 remaining objects, only 
one object (their ID#70) shows a sign of emission-line ex- 
cess (green open circles in Fig. 13). However, the object 
is just below our HAE selection criteria (NB< 22mag and 
E > 2) and therefore not present among our HAE candi- 
dates. 

We comment on the failure of the HAE counterpart 
detection for their Lya emitters in detail. We can calcu- 
late the expected Ha fiux from their cataloged equivalent 
width, assuming the Lya/Ha ratio of 8.7 for the Case B 
recombination (Osterbrock 1989). We find that the ex- 
pected Ha flux are mostly several xlO""'^'' erg cm~^, 
close to our detection limit. However, this should be seen 
as a lower limit, as the Lya/Ha can be much smaller 
than 8.7 due to the more effective dust absorption by the 
resonant scattering nature of the Lya line (Kurk et al. 
2004a; Hayes et al. 2010). Thus at least some, if not all, 
of their 9 LAEs could be detected. We note that the cat- 
aloged if -band magnitudes of Lya emitters are all quite 
faint (5 are upper limits with I^ab > 24.35, and only two 
are within our magnitude limit of i^AB > 23.5). The non- 
detection by our K image could be reasonable. However 
they arc also totally absent even in our NB2288 image. 

Kurk et al. (2004a) also did not detect Ha emission 
for their sample of Lya emitters in the PKS 1138 proto- 
cluster. For two of their three brightest Lya emitters, for 
which Ha emission could be detectable, faint emission was 
seen, but not strong enough to pass the HAE selection cri- 
teria. The range of the upper limit of the Lya/Ha ratio 
for their Ha emitters arc 0.03— 1.06, significantly lower 
than the Case-B assumption (see also. Hayes et al. 2010). 

There are two possibilities to explain this. One possibil- 
ity is that the redshift of the Lya emitters are out of our 
NB2288 narrowband window. The fiher used by KC97 
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has a central wavelength of 4229A with a filter width of 
70A. The Lya emission that can enter into the window 
should be in the range 2.45 < z < 2.50, and the corre- 
sponding Ha line will come to 2.261 fim~ 2.294 /im. As 
is shown in Fig.l, almost half of the shorter wavelength 
range (2.261 //m — 2.28 ^m) is out of our survey window. 
However, if the Lya emitters are all such objects, they 
would not be associated with, but slightly foreground to 
4C 23.56. 

Another explanation is a large error in the Lya EW 
estimate, or spurious detections. Due to the lack of B- 
band data which samples the continuum just redwards of 
Lya, they had to use a crude UV continuum estimate. 
If the EW of their Lya emitters are much smaller than 
their quoted values, our non-detection of Ha counterparts 
could be understood. An independent Lya emitter search 
is necessary to check the reliability of their selection of 
NB excess objects. We are in the progress of carrying out 
such a Lya imaging survey for this field. Its results will 
be published elsewhere (Tanaka I. et al., in preparation). 

A. 2. 2. Objects in Kajisawa et al. (2006b) 

Kajisawa et al. (2006b) took the JHK images of the 
4C 23.56 field and found a significant excess of ^^JHK' 
objects (see also, Kodama et al. 2007). We checked their 
Ha excess in our catalog. Though we could not identify 
their b- JHK objects, several Distant Red Galaxies (DRGs) 
were identified. The color excess of their DRGs is shown 
in Fig. 13 (blue filled circles). Excluding 4C 23.56, two 
objects are identified in our HAE sample. Another two 
objects are identified as KC97's EROs. One of them shows 
emission- line excess (#70 in KC97). Two other objects do 
not show detectable emission-line excess. 

A. 2.3. An Object in Stockton et al. (2004) 

Stockton et al. (2004) reported an object in 4C 23.56 
having very red color and a disk-like morphology. The 
SED fitting result indicates the strong preference for the 
single-burst model with fairly old age of ~ 2Gyr at the 
redshift of the radio galaxy. The result suggests that the 
galaxy should be inactive Stockton & McGrath (2007), 
and no Ha emission is expected. 

Their object is actually identified as an ERO in KC97 
which is a DRG in Kajisawa et al. (2006b). In Fig. 13 we 
show the color excess of the object (magenta open trian- 
gle). As expected, the object does not show signs of excess 
emission from Ha. 
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Fig. 1. (Top) The simulated shift of the transmission curve of the NB2289 filter. The horizontal axis is the wavelengths in /im and 
the vertical axis has arbitrary units. The curves with orange (dot-dot-dashed), red (thick solid), magenta (thick solid), and blue 
(thin solid) colors each indicate the transmission curves at x=256 (close to the innermost FOV coverage), 1024 (detector center), 
1536, and 2048 (outer edge of the FOV) in channel-2 raw-image coordinates (the bottom figure, and corresponds to the right-half 
of the MOIRCS FOV). At each x position, the shift of the transmission curve along with {/-coordinates is also shown: the curve at 
the longest wavelength is at the detector center (j/=1024), and shifts to shorter wavelengths as the y coordinates goes to the edge 
of the detector (y=512/1536:middle and l/2048:bluest). The red arrow indicates the expected position of the Ha emission, and 
the vertical green lines (dots) indicate the range of the shift for objects at ±500 km relative to the radio galaxy. The peak 
amplitude of the transmission curves at each x coordinate are artificially shifted for clarity. See discussion for §2.1.1 for details. 
(Bottom) The schematic figure of the MOIRCS field of view for explanation of the top figure. The full 4' X 7' FOV is divided into 
two channels. Note that the minimum x coordinate of the raw data area is around x = 270 for channel 2 (x < 270 is blocked by the 
mirror shadow). The curves on the channel-2 side show the iso-redshift contours (each number indicates the redshift) we can detect 
Ha at the central wavelength of the filter window. The NB2288 filter can cover the redshift range of ±0.017. The pointing center is 
around (270,1024). 
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Fig. 2. The ii's— NB versus NB color-magnitude diagram for our cataloged objects (plus with the label "ALL"). The blue asterisks 
(labeled as "EM" ) indicate our selected emission-line objects based on the criteria EWo > 50 A, the narrowband aperture magnitude 
of less than 22 mag, and S > 2, which is used for analysis in the paper (one spurious object in the sample is removed here). The 
brightest emitter is 4C 23.56 which is shown as a large asterisk. Emitters with red open boxes are the objects with MOIRCS spectra 
(except 4C 23.56). The tilted 4-dashed line indicates the 3-sigma limiting magnitude of Kb = 23.5. The insetted figure shows the 
transmission curve of the filters used in this study. 
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Fig. 3. MOIRCS MOS spectra of four Ho emitters. The top spectra (HAE#491) is 4C 23.56. The barcodc-likc pattern at the 
bottom is the OH airgrow spectra. The numbers at the bottom are wavelengths in angstrom. 



20 



Tanaka et al. 




22000 22500 23000 23500 22000 22500 23000 23500 
Wavelength (Angstrom) Wavelength (Angstrom) 




22000 22500 23000 23500 22000 22500 23000 23500 
Wavelength (Angstrom) Wavelength (Angstrom) 



Fig. 4. MOIRCS MOS spectra of four Ha emitters. The top left spectra (HAE#491) is 4C 23.56. The OH airglow lines (red solid 
line) and the atmospheric transmission curve (blue dotted) arc also plotted at the bottom panel of each spectra. Note that the 
spectra are not calibrated. 
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Fig. 5. (a: left) The B — J versus J — /^s two-color diagram for objects in the INGRID subfield. Five red boxes with the label 
"HAE(OKAY)" indicate the HAE candidates for which we can measure colors. Small plus symbols indicate other NB-selected 
objects in the area. The overplottcd model tracks from blue to red indicate the "pure disk" model with continuous SFR (bluest), 
"50% disk + 50% bulge", "20% disk + 80% bulge", "pure bulge" (reddest, and mostly out of the figure) models, respectively. The 
cyan-colored region on each model track indicates the range 2.0 < z < 3.8, and the farthest point in each cyan-colored track is around 
z = 2.5 to 3. See text for more details, (b: right) The same figure for galaxies in GOODS-North field (see Appendix for details). 
The data is from the MOIRCS Deep Survey (MODS) Catalog with Ks < 23 mag (Kajisawa et al. 2009; Kajisawa et al. 2010). In the 
figure, red filled circles indicate objects with spectroscopic redshifts of 2.3 < z < 2.6, similar redshifts to our Ha emitters at z ~ 2.5. 
Blue triangles indicate objects at 0.2 < z < 0.25 (open) as well as 0.75 < z < 0.8 (filled), each of which corresponds to the redshift 
range for our foreground Paschen series (Paa and Pa/3, respectively) interlopers, and the filled boxes (skyblue) indicate objects at 
3.3 < z < 3.8 to which the background [O III] emitter contaminant for our data would belong. Inset histogram shows the fraction 
of objects with known spectroscopic redshift that satisfy the BJK criteria we used here as the object selection at z > 2. Note that 
there is a very limited number of spectroscopically confirmed galaxies at z > 2.7, as indicated by the large N^/^ errorbars. 
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Fig. 6. The distribution of the HAEs on the NB2288 image. The FOV covers 7' X 4' on the slty, with north and cast to the top 
and left, respectively. The red open circles indicate the candidate HAEs by our selection. The three small open boxes indicate the 
Distant Red Galaxies from Kajisawaet al. (2006a) that show emission-line excess (i^s— NB> 1) as seen in Fig. 13 (also see Appendix). 
Two of them are also identified as HAEs. The remaining one is the only NB-cxcess object among the ERO sample by KC97. Two 
small open circles (green) indicate the candidate emission-line objects that satisfy the criteria of S > 2 but are fainter than NB=22.0. 
Three red open circle with the letter "S" indicate the emitters with the MOIRCS spectra (§3.1.1). 4C 23.56 is shown as a large 
diamond. The area enclosed by the white box indicates the same area shown in Fig. 7. The yellow rectangular region indicates the 
INGRID subfield (see §3.1). 
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Fig. 7. Zoom-up views of the "Eastern clump" region in both the NB2288 (right) and the Ks (loft) images. The objects at the 
center of each yellow circle represent the HAEs. The Icr, 4(t, 7a", lOcr, &:13(t contours based on the rms noise of the Ks image before 
smoothing is applied to the Gaussian-smoothed Kb images, and is also overlaid on the NB2288 image. A clear extended or offsettcd 
light distribution of HAEs in the NB2288 image is seen for (at least) the three objects on the right. Interestingly, the directions 
of the offsetted light distributions for these three objects are roughly the same (to the North- West). The size of the area is about 
26" X 26" (~ 0.5 X 0.5 Mpc^ comoving). 
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Fig. 8. The difference of the object position between the Ks and the NB2288 images. Objects that have similar magnitudes to our 
HAEs are used (dots). The 1-, 2-, and 3-(t scatter of the distribution are shown as dashed circles in the figure from the inner to 
outer, respectively. The result for the HAEs is shown by the filled squares, with the 4 HAEs in eastern clump in Fig. 7 overlaid by 
asterisks. The offset is shown in X and Y pixel coordinates, where the MOIRCS pixel scale is 0"ll7/pixel (corresponding physical 
scale is 0.96 kpc/pixel at z = 2.48). 
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Fig. 9. The density contour of the faint Spitzer MIPS 24/xm sources is overlaid on the MIPS 24/im image. The area with an 
exposure coverage less than 1/3 of the total is not shown here. The brightest source in the center is 4C 23.56. The red rectangular 
region indicates the MOIRCS field of view (North is up: see Fig. 6), and the red open circles indicate the HAEs. The distribution 
of the faint MIPS sources (90 < -F24/im [/^Jy] < 275) is smoothed with a Gaussian kernel with o" = 30", and shown as contours. The 
contour interval is (2i + l)-times (i = 0, 1, 3, 6, 11) the square root of the average MIPS source density (from the "outer region" : see 
text) above average. 
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Fig. 10. The Ks versus stellar mass relation for objeets at 2.3 < z < 2.6 from the MOIRCS Deep Survey (MODS) eatalog used in 
Kajisawa et al. (2009). The eatalog ineludes the eompilation of the spectroseopic redshift data published, and are shown as crosses 
in the figure. The stellar mass is based on the SED-fitting by the 12-band photometry from UV to infrared, with the SED model 
from Bruzual & Chariot (2003). Details of the model fitting procedure can be found in Kajisawa et al. (2009). 
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Fig. 11. (a: left) The observed star-formation rate (SFR) of the HAE sample as a function of stellar mass. The plus marks indicate 
the Ho-based SFR before the extinction correction, while red filled box connected to each plus mark is its extinction-corrected SFR 
values (see text). The blue filled boxes are the SFR based on the SpitzerMlPS 24 fim photometry. The dotted line (blue) indicates 
the mass-SFR relation at z ~ 2 from Daddi et al. (2007). The two arrows indicate the objects with the upper-(lower-)limit for their 
mass (SFR before extinction correction) due to their faintness. The most massive object is 4C 23.56. Its MlPS-based SFR is not 
shown (due to an AGN contribution and the model limit), (b: right) The comparison of the Ho-based specific SFR (SSFR) of 
our HAE sample and of some field Ha-emitting objects as a function of stellar mass. The marks with dust extinction correction 
(red filled box) and raw values (plus) are connected by vertical lines. The field samples are from Forster-Schreiber et al. (2009) and 
Onodera et al. (2010). Note that Forster-Schreiber et al. (2009) shows the results using two types of the extinction correction; thus 
they are shown as connected vertical lines in the figure. Our SSFR estimate should be treated as the crude values due to the large 
uncertainty in both the extinction correction and the mass estimate. Note that 4C 23.56 is the point at A/ > 10^^ Mq, and should 
be ignored because of strong AGN contamination in both the K-hand and the Ha emission. 
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Fig. 12. Comparison of the Kb magnitude and L(Ha) for Ha emitters from two field surveys from Forster-Schreiber et al. (2009) 
and Onodera et al. (2010) (asterisk and cross), PKS 1138-262 field from Kurk et al. (2004a) (box), and our study (plus). Only 
objects at z > 2 are shown for field survey data. For the PKS 1138 field, wc plot the data with the significance level over 3 and 
EWo > 50A for better match to our selection criteria. Two arrows at brighter than K = 20 represents the PKS 1138 and 4C 23.56 
radio galaxies which are out of the figure due to their bright magnitudes. 
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Fig. 13. Color excess of objects studied by others in the 4C 23.56 field. Only objects identified in our catalog are shown here. 
The grey plus marks are all objects in the field of view, and the asterisks are the HAEs. A red filled box and 4 green open circles 
each indicate a Lyo emitter and the EROs from KC97 that are in our catalog. Blue filled circles indicate the DRG sample from 
Kajisawa et al. (2006b), and a magenta open triangle is the object studied by Stockton et al. (2004). Ignoring the radio galaxy, 
three DRGs/EROs show the emission-line excess [Kg — NB > 1), and they are also plotted in Fig. 6. 
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Table 1. Observational properties of candidate Ha emitters. 



HAE ID 


RA(J2000) 


DEC(J2000) 


Kj 


EWo^ 


F(Ha)-^ 


F(24 ^m)-* 


redshift'^ 


153 


21:07:16.96 


23:30:30.5 


21.82 ±0.07 


327 


3.64 ±0.77 


78.9 


2.4879 


354 


21:07:21.03 


23:31:13.9 


21.23 ±0.05 


105 


2.47 ±0.76 


134.9 


2.4865 


356 


21:07:15.70 


23:31:11.9 


22.78 ±0.12 


282 


1.35 ±0.50 


125.9 




393^^ 


21:07:21.48 


23:31:19.5 


21.72 ±0.06 


51 


0.81 ±0.55 


42.4 




431 


21:07:15.93 


23:31:27.5 


22.17±0.08 


221 


1.96 ±0.62 


156.6 




479^ 


21:07:21.63 


23:31:41.2 


22.37 ±0.08 


956 


4.25 ±0.72 


131.8 




491*^ 


21:07:14.80 


23:31:45.1 


19.49 ±0.02 


122 


13.95 


4630 


2.4852 


500^ 


21:07:21.84 


23:31:41.8 


> 23.5 


> 780 


>0.8 






51110 


21:07:22.58 


23:31:43.4 


22.98 ±0.12 


300 


1.18±0.46 


156.5 




526" 


21:07:21.78 


23:31:49.6 


20.76 ±0.04 


95 


3.45 ±0.94 


119.2 


2.4872 


543 


21:07:17.61 


23:31:50.5 


> 23.5 


> 505 


>0.8 







1. the SExtractor MAG_AUTO output in AB mag. 

2. rest-frame equivalent width based on the (ii:s-NB2288) values by MAG_AUTO output. 

3. continuum subtracted narrowband flux in 10"^^ erg sec~i cm~^. 

4. Spitzer MIPS 24 ^m flux in ^Jy from the DAOPHOT photometry. 

5. see §3.1. 

6. a possible low redshift interloper based on the BJK criterion (§3.1). 
8. 4C 23.56. Its magnitude is uncertain due to stray light (see §2.1). 
7,9,10,11. galaxies in the "Eastern Clump" (Fig. 7; sec §3.3). 



Table 2. The estimated basic properties of the Hq emitters. 



HAE ID 


log(M/M0) 






gpj^corr. 


l0g(i24 Mm/L©) 


SFR(24 ^iuif 


153 


11.124±0.235 


14.4±3.0 


67.13 


187.96 


12.533 


410 


354 


11.430 ±0.225 


9.75 ±3.02 


45.57 


355.45 


12.792 


743 


356 


10.635 ±0.260 


5.34 ±1.96 


24.94 


122.21 


12.758 


688 


393 


11.178 ±0.230 


3.19 ±2.16 


14.92 


14.92 


12.237 


207 


431 


10.948 ±0.243 


7.74 ±2.45 


36.17 


177.23 


12.864 


877 


479 


10.844 ±0.240 


16.75 ±2.82 


78.30 


383.67 


12.781 


724 


491" 


< 12.320^ 


55.0 ±6.8 


257.15 


1429.75^ 


N/A^ 


N/A8 


500 


< 10.164 


>4.63 


> 24.85 








511 


10.531 ±0.262 


4.65 ±1.8 


21.72 


106.43 


13.091 


876 


526 


11.667 ±0.220 


13.6±3.7 


63.67 


311.98 


12.869 


647 


543 


< 10.164 


>4.13 


> 19.28 









1. Extinction- uncorrected values in unit of 10 erg s ^. 

2. in M© yr~i. See text for detail of the extinction correction. 

3. Spitzer MIPS 24 /xm-based SFR in Mq yi'^ . 

4,5. 4C 23.56. Its mass estimate is unreliable due to the uncertainty in its Kg magnitude (see §2.1). If we use the 

value by KC97, the estimated stellar mass is 11.67. 

6. Affected by an AGN contribution. The SFR value is quite uncertain. 

7,8. Affected by an AGN contribution. The L24 ^ra conversion from the observed MIPS flux is impossible due to the 
limitation of the Chary and Elbaz (2001) model. 



